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Abstract 

In this note we consider a system with a large angular momentum i 
whose state we can store using some log 2 (i') qubits. The problem then 
is how to carry out spatial rotations of the system in this representation. 
In other words we are looking at a unitary representation of SU(2) with 
dimension 21 + 1 and want to implement these transformations with re- 
sources polynomial in log(^). We only give a sketch of our solution which 
involves "storing" discretised spherical harmonic functions Yi m (®,<f>) in 
a quantum register. Also there are some technical gaps in the construc- 
tion, but they are based on plausible assumptions. Our approach is rather 
cumbersome and we hope somebody will find a nicer solution. For a nice, 
elementary explanation of what we are trying to do (not involving physics 
or representation theory) see section 14.6.21 

1 Introduction 

Consider a system (say a particle) with total angular momentum given by £, 
thus the absolute value squared of the angular momentum is £{£ +1). The 
21 + 1 dimensional state space (Hilbert space) of this system is spanned by 
the eigenstates of the z-component of the angular momentum with eigenvalues 
r7i = — i ... + £. (For the time being we only consider integer £, thus e.g. orbital 
angular momentum in quantum mechanics.) We can now easily represent a 
(pure) state of this system in a quantum register with some log 2 (^) qubits, 
whereby we choose the m-eigenstates to correspond to the computational basis 
states: \m), m = —£... + £. 

A rotation about the z-axis is easily implemented as |m) — > e lm ^\m). If we 
knew how to implement some (actually almost any) other spatial (thus SO(3)) 
rotation we could combine it with z-rotations to obtain any rotation. But doing 
such a rotation has proved quite tricky and is the main point of this note. 



In more mathematical terms we are simply trying to efficiently implement the 
unitary transformations corresponding to a high dimensional representation of 
SU(2) (or SO(3)) on a quantum computer. 

The approach we present here is rather cumbersome but it seems to do the 
job of efficiently implementing any rotation with arbitrarily small (exponentially 
small) error. Still there are some gaps, but we believe these are just (tedious) 
technicalities. Also, as mentioned before, we hope that someone finds a much 
nicer (possibly exact) solution, although this has eluded us. 

2 The approach 

The idea is that smooth functions in 3 dimensions (thus functions from R 3 — > C) 
form an SO (3) representation under spatial rotations. In particular the spher- 
ical harmonics Yf m (0,</>) for a fixed i transform under spatial rotations like 
the angular momentum eigenstates thus they form an irreducible rep- 

resentation of SO(3) (and thus also of SU(2)). (Thus a rotated spherical har- 
monic can be written as a linear combination of Yf m 's with the same £.) Thus 
in a way we are simply looking at a wavefunction of a particle in 3 dimen- 
sions which has the right (orbital) angular momentum. (Because for the an- 
gular momentum the radial dependence doesn't matter, we can actually look 
at wavefunctions which depend only on 9 and </>.) For a nice picture of a 
random superposition (drawn from a Gaussian distribution) of spherical har- 
monics with £ = 20 look at http://www.zalka.itp.unibe.ch/ylm.jpg or 
http : //www . iqc . ca/~zalka/LANL/ graph/ylm . jpg . 

The idea now is to represent ( "store" ) such functions on the quantum com- 
puter in a discretised form using a sufficiently fine 3-dimcnsional cubic lattice. 
Thus we have 3 quantum registers \x, y, z) which we prepare in a superposition 
with amplitudes equal to the function value at the point x, y, z (apart from an 
overall normalisation) . Still, we may ultimately like to apply the (representation 
of the) rotations to the compact encoding in a single register |m). For this we 
first have to translate this encoding into the more complicated one represent- 
ing a 3-dimensional function. Thus to implement a rotation we have to do the 
following 3 (unitary) steps: 

• translate: \m) — > \Yg m ) 

• rotate: spatially rotate \Yi m ) 

• translate back: \Ye m ) — > |m) 

Where everything should just as well work for superpositions over different 
m's. The state \Yi m ) stands for a suitably discretised version of the function 
Y( m (Q,(j>) on a 3-dimension grid in the 3 registers \x,y,z). The discretisation 
has to be fine enough to faithfully represent the smooth function, but this can 
be done with only 0(\og(£)) qubits. We also have to choose some r-dependence 
and I propose to simply consider the spherical harmonics on a thin spherical 
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"shell". Thus really we have the overall function Y^ m (G, <j>) ■ R(r) where R(r) is 
only non-zero in a narrow region around some ro where e.g. we choose it to be 
constant. 

One way to do the "translation" is to first think how we could prepare the 
state \Yi m ) for a fixed m. If we knew that (and I will sketch how that could be 
done), we can implement the unitary transformation \m) — > \m,Y( m ). To see 
how we then can get rid of the m, consider the converse: given |Yf TO ), how can 
we find to? For this, note that a spatial rotation by <f> about the z-axis of Yg m 
gives an overall phase of e lm ^. Thus, if we know how to spatially rotate |Y« m ), 
we can do "phase estimation" (see Kitaev, Mosca) to do |Yf m ) — > \Y^ m ,m). 
Combining these techniques approximately gives a unitary transformation doing 
|m) -» \Y lm ). 

How can we do the spatial rotation of \Ye m ) ? Of course rotations by 90° 
about the principal axes are easy as the points on the rotated grid (a cubic 
lattice) match those of the original one, e.g. we can do \x,y,z) — > \x,z,—y) 
reversibly and thus unitarily. For rotations about some generic angle we have to 
find a 1-to-l mapping between the points of the rotated and the original lattice. 
Small errors, thus mapping a point to another point a few lattice spacings away 
won't matter, as we can make the lattice spacing exponentially small and thus 
the effective error will be small. Unfortunately simply mapping a point to the 
nearest point of the other lattice isn't 1-to-l and thus not unitary. Thus we 
have to think of something different (see below). 



3 (some) techniques 
3.1 rotation on a lattice 

What we need is a 1-to-l mapping between points of two cubic lattices which are 
rotated with respect to each other. As it is sufficient to do rotations about one of 
the principal axes, the problem really reduces to two 2-dimensional (square) lat- 
tices with the same origin. The idea is to do several "shearings" which together 
give a rotation. A "shearing" is a linear map of the form x — > x + a ■ y, y — ► y, 
thus a mapping where all lines parallel to one of the axes are mapped into 
themselves. Doing 3 such mappings: 




we can e.g. determine b in terms of a such that overall we get a rotation 
(this is an easy exercise) . The individual "shearings" can rather easily be done, 
as they simply map lines into themselves with a translation. We can then easily 
compute a map (and its inverse) that e.g. maps each point on a line to the next 
one of the other lattice. 

But the lines are finite (as the lattice is finite). Then, say, we use periodic 
boundary conditions, thus our mapping is cyclic. At any rate, the outer parts of 
the whole lattice will not be mapped according to a rotation, this only happens 
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(approximately) for a region around the centre of the lattice. But this is no 
problem, as we simply make sure that our wavefunction is only non-zero in this 
inner region. 

3.2 preparing spherical harmonic states |Yg m ) 

Thus we want to prepare a state (of a quantum register) whose amplitudes 
are taken from a smooth, "well known", natural function. Let's first consider 
the 1-dimensional case. Thus we have a "natural" function f(x) and want to 
prepare the state ^ x f(x)\x) (assuming this is normalised) where, say, x = 
... 2™ — 1 . In this state the probability of finding the highest order bit in state 

|0) is J2x=o 1 \f( x )\ 2 > which we can achieve by suitably rotating (say with an 
SO(2) rotation) a qubit initially in state |0). To get the right probabilities 
for the TWO high order bits we have to suitably rotate the second highest 
qubit (initially in state |0}). First we have to determine the sums over the 4 

subintervals Extf ^ 1 l/<»| 2 with a = °) 2 "^ 2 ; 2-2"~ 2 and 3-2™" 2 . Depending 
on the values of these 4 sums and conditioned on whether the highest bit is or 
1, we have to rotate the second highest qubit. And so on. Thus at one point we 
have to rotate qubit number i by an angle which depends on the values of all 
higher bits (and on the function). The correct rotation angle will be given by 

the ratio of two sums of the form Y^x=k\? * l/( x )| 2 with k — fco an d k = fco + 1. 

Thus before carrying out the controlled SO(2) rotation of qubit number i, 
we have to compute the value of two of these sums (which for smooth functions 
can be approximated by integrals). This computation has to be carried out 
"in quantum parallelism" . Thus the value of the higher order bits determines 
which sums we have to compute. This (approximative) computation of integrals 
is the most cumbersome part of our construction. Still, for "natural" , "well 
known" functions like the spherical harmonics, we can expect to find an efficient 
algorithm for this. In the end (if we use conditional SO (2) rotations), we get a 
state that has the right amplitudes, possibly up to phase. But a "rephasing" , 
again according to a well known, natural function is rather easy to carry out 
efficiently. This whole technique is e.g. described in section 3.2 (probably 
more carefully than here...) and has also been described by other authors. 

Clearly the technique also generalises to functions of several variables (lead- 
ing to a state of several quantum registers). Instead of integrals over "binary 
intervals" we will now have to compute integrals over "binary hyper-intervals" . 
In our case we have the function Yg m (®, <fi) ■ R(r) on a 3-dimensional lattice. In 
a first step we could get the right probability distribution over z-values. Thus 
in the registers \x,y, z) we would have x = y = and thus a distribution con- 
centrated on the z-axis. For this first step we would only have to consider the 
generalised Legendre functions Pe m (z — cos(O)) and integrals over it. (Remem- 
ber that Yfo„(0, 4>) — P^ m (cos(0)) • e lm ^.) From there all we would have to do 
is (thinking of one z = const, plane) to "spread out" the x = y = point into 
a ring and do a rephasing according to e 4m< ^. (It's a ring because we chose a 
narrow, rectangular radial function R(r).) 
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The main technical gap in our construction is to show how integrals over 
the generalised Legendre functions Pe m (z) can be computed (approximatively) 
efficiently. I believe it is clear that this is possible! 

4 various remarks 

4.1 How to measure I? ?! 

We have shown how m (z-component of the orbital angular momentum) could 
be measured for a wavefunction on a lattice. But I don't know how to measure £, 
thus the absolute value squared of the angular momentum. L 2 = L 2 +L 2 +Ll = 
Ae,0- But even though we know how to measure L 2 , L 2 , L\ (and also A Xj y ; z) 

this doesn't allow us to measure L 2 , because these terms don't commute. If we 
knew how to measure L 2 , this would allow us to initialise the \x,y,z) register 
in a specific irreducible representation of SO (3), although maybe not one of our 
choosing. Still, for many applications this would be enough and would relieve 
us from doing the cumbersome \m) — > \Yi m ) transformation. Failing other 
ways, one way of measuring L 2 might be to do a "full translation" \£, m) — ► 
\Yi m ) (involving estimating I from |5^, m )) but then that would be again very 
cumbersome. 

4.2 Doing without the "compact representation" ? 

As mentioned above, sometimes we may not need to apply the transformation 
to the "compact" \m) representation, just applying it to the 3D lattice repre- 
sentation might be good enough. (E.g. if all we need is some kind of "phase 
kickback", it doesn't matter what representation we have.) Still, we may want 
to initialise things if not in a specific state, then at least in a specific irreducible 
representation. But, as mentioned above, we haven't found an easy way of doing 
this. 

4.3 Doing things on a 2D ©, (ft lattice. 

There are alternatives to representing the spherical harmonics Ye m (Q,<j>) on a 
3-dimensional cubic lattice. E.g. we could simply use a 2-dimensional cos(0), <f> 
plane (a square lattice). The problem is that the spatial rotations then get 
more complicated. E.g. a rotation by 45° about the x-axis (which would be 
sufficient) would correspond to some smooth (but not linear) mapping of this 
plane into itself. Again we could try to work out how to make a 1-to-l mapping 
to close points (the "smooth mapping" is area preserving, thus the density of 
points remains unchanged), possibly by piecewise linearising the smooth map. 

An advantage of using the cos(0), <f> plane would be that we wouldn't have 
the problem of uncomputing m: |Y^ TO ,m) — ► \Yi m ). We could first construct 
Pt,m(z)'- \ m ) ~ > \Pe,m(z),m) and then simply Fourier transform m to get the 
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correct ^-dependence. Possibly similar techniques (avoiding the "phase estima- 
tion") are also possible for the 3D lattice. 

4.4 "Changing lattices": useful for quantum simulations. 

So we represent a (smooth) wavefunction on a lattice in some coordinate system. 
Sometimes we may want to change the lattice, e.g. by rotating it, as we did here, 
or e.g. by changing from a Cartesian lattice to a spherical coordinate one (with 
same density of lattice points). This is similar to a problem in image processing 
when we want to change the grid of pixels (we may call this "rescanning"). But 
here, in order to do things unitarily, we have to find a 1-to-l mapping between 
close points. (Unless we want to do thing non-classically, thus mapping a point 
on one lattice to a superposition of nearby points of the other lattice.) E.g. if we 
want to measure the z-component of the angular momentum, a polar coordinate 
lattice would be nicer. (This would e.g. give a nicer way of estimating m from 

\Yt,m)-) 

4.5 Why don't "efficient quantum simulations" solve our 
problem ? 

E.g. Seth Lloyd (Science, August 1996) (but also [T|) has shown that physical 
quantum systems can efficiently be simulated on a quantum computer. One 
may wonder whether that would not already include what we are doing here, 
namely "efficiently simulating large angular momenta" . The point is that it has 
been shown how to do simulations that are polynomial (actually linear) in the 
physical size of a quantum system, thus space, time, energy, precision, etc. But 
although we can store the state of an angular momentum I using only 0(\og(£)) 
qubits, the physical size (here space, energy) of such a system is necessarily 
large, namely polynomial in £. Thus what we are trying to do is somewhat 
more ambitious than an "efficient quantum simulation" . 

Apart from angular momentum one may think of simulating other "few- 
particle" system with high dimension, e.g. two interacting harmonic oscillators 
(truncated at some high energy). 

4.6 "Adding" angular momenta 

"Angular momentum addition" is a standard subject in quantum mechanics (in- 
volving Glebsch-Gordon coefficients etc.). Two systems (thus formally a tensor 
product) with angular momentum quantum numbers I and £' carry irreducible 
representations I — £' to £ + £' (assuming £ > £'). Thus in a way if we have 
representations £ and £' we have in particular also £ + £'. It is also not difficult 
to find a "translation" from a compact encoding of the large representation (on 
a QC) to the two smaller ones: \M) — * c ■ J2 m l TO ' m' = M — m). Thus if we 
know how to simulate I and £' , then we can also do I + £'. 

This is why our problem would already be solved even if we could only "pro- 
duce" representations of a random size, or only dimensions equal to a power of 
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2. Also the difference between integer and half- integer (spin) angular momenta 
is no problem, as we can always "add" a single qubit (a spin 1/2 representation). 

4.6.1 SU(2) or SO(3) representations? 

For integer t we have a representation of both SU(2) and SO (3), while for half- 
integer angular momentum (thus spin) it is only a representation of SU(2). Re- 
member that there is a natural 2-to-l homomorphism from SU(2) to SO(3), thus 
they are "almost the same". With the spherical harmonics Y^ m (0,</>) we have 
all SO(3) representations (all odd-dimensional SU(2) representations). Proba- 
bly there is a similar representation (in terms of smooth functions) that directly 
gives all SU(2) representations. But spherical harmonics are more familiar to 
me. 

4.6.2 Symmetric subspace of TV qubits. 

There is a nice way to understand what we are trying to do here for people not 
familiar with angular momentum (physics) or representations (of Lie groups). 
Consider N qubits. Now consider the space of (pure) states invariant under any 
permutations of the qubits. This "symmetric" subspace is clearly spanned by the 
"uniform amplitude superpositions" of all binary strings with equal Hamming 
weight. Thus one basis state in this symmetric subspace is the superposition 
(with equal amplitudes) of all strings with h l's (and N — h O's). So for each 
Hamming weight we have such a basis state, thus the dimension of the symmetric 
subspace is (only) N + 1. 

Now imagine we apply the same transformation to each of the N qubits (thus 
overall U® N with U G SU(2)). Note that if we start with a state in the "sym- 
metric subspace" we remain in it, thus this is an invariant subspace. What we 
are trying to do in this note is to "efficiently" encode the symmetric subspace in 
some log 2 (N) qubits and carry out the transformations corresponding to doing 
the same thing to all N original qubits. We already know how to do transforma- 
tions corresponding to diagonal SU(2) matrices. To be able to do everything, 
it would clearly be sufficient to e.g. do the transformation corresponding to 
doing the Hadamard transformation on each of the N qubits (thus H® N ). Note 
that this "hypcr-Hadamard" can easily be written down by thinking in terms 
of the symmetric subspace (use H® N \x) = 2~ N / 2 J2 y (— ^) x ' v \y) and a bit of 
combinatorics) . 

What we are doing with the "symmetric subspace" of N qubits simply cor- 
responds to "adding" N spin 1/2's in the above sense. 

4.7 Representations of Lie groups other than SU(2). 

It seems that once we know how to "implement" SU(2) representations, we can 
use this to do any Lie group representations (e.g. SU(n), Sp(n), Es, etc.). This 
is because in some sense these groups (the "continuous groups") can be built 
up from SU(2). 
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4.8 Using this to simulate the chaotic "kicked top". 

This investigation was motivated by the idea to simulate the (quantum-) chaotic 
"kicked top" on a quantum computer (related to me by Raymond Laflamme and 
David Poulin, Waterloo). There one considers a "large spin" j and repeatedly 
applies to it a rotation about the y-axis and a diagonal (unitary) transformation 
of the form \m) — > e lc ' m \m). The latter transformation is easy to implement 
(as is any rephasing e 1 ^ 171 ^ if we can compute /(m), see e.g. eq.(7)). Thus 
all that is missing is the y-rotation, which is an SU(2) representation. 

4.8.1 Doing it on an NMR quantum computer. 

An interesting simulation of this "kicked top" could also be carried out on 
an "ensemble" quantum computer, thus a quantum computer where only a 
constant number of qubits is properly initialised in state |0) while the rest are 
in a maximally mixed state. (This model corresponds to NMR.) The problem 
is that in our construction, to carry out computations, we need a lot of properly 
initialised ancilla qubits. Thus it would be nice if a future "nice and natural" 
implementation of SU(2) representations could do it without using any such 
ancilla qubits (like the Quantum Fast Fourier Transform QFFT), or then only 
with a constant number of them. 

4.9 Any connection to the Fourier Transformation? 

The SU(2) representation transformations form a natural class of unitary trans- 
formations, a bit like the Fourier transformations (say for cyclic groups). In 
particular we can look at the representation of the Hadamard transformation 
(a kind of "hyper-Hadamard" ) and ask whether it has any connection to the 
Fourier transformation. We haven't found any such connection, except maybe 
in the limit of large dimensions. 
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